Analytic solution of the separability criterion for continuous variable systems 
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By using the algebraic separability criterion of R. Simon, we present an explicit determination 
of squeezing parameters for which the P-representation condition saturates the Sp{2, R) Cgi Sp{2, R) 
invariant separability condition for continuous variable two-party Gaussian systems. We thus give for 
the first time the explicit analytic formulas of squeezing parameters which establish the equivalence 
of the separability condition with the P-representation condition. The implications of our algebraic 
analysis on some of the past related works are discussed. 

PACS numbers: 



I. INTRODUCTION 



The entanglement is a basic notion in quantum me- 
chanics, but the quantitative criterion of entanglement 
is known only for a simple system such as a two-spin 
system [H, In view of this fact, it is remarkable 
that a proof of the necessary and sufficient separabil- 
ity criterion for continuous variable two-party Gaussian 
systems has been given by R. Simon [31 on the basis 
of generalized Peres- Horodecki criterion l|, He also 
gave an algebraic criterion for separability (i.e., non- 
entanglement) 0, though it was not used in his proof. 
This algebraic criterion is considered to be fundamental, 
but its explicit analysis has not been performed so far. 
See a review [3] on the present status of the quantum 
separability problem of Gaussian systems. 



We here present an explicit determination of squeez- 
ing parameters for which the P-representation condition 
saturates the Sp{2, R) (g) Sp{2, R) invariant separabil- 
ity condition by explicitly solving the algebraic condi- 
tion of Simon. We thus give for the first time the ex- 
plicit formulas of squeezing parameters, which establish 
the equivalence of the separability condition with the P- 
representation condition, in terms of the parameters of 
the standard form of the covariance matrix (or second 
moments) for Gaussian systems. These explicit analytic 
solutions should be useful in the quantitative theoretical 
and experimental analyses of entanglement such as in ^ . 



We also show that our analytic solutions of squeezing 
parameters ri and r2 do not satisfy in general the equa- 
tion f{r^) — which appears in another formulation of 
the separability criterion for two-party Gaussian systems 
In this sense our scheme is quantitatively different 
from the scheme in ^6j. ft is however shown that our 
exlicit analytic solutions allow us to construct a concrete 
proof of the separability criterion on the lines of Q. 



II. ANALYTIC SOLUTIONS 

We start with the 4x4 correlation matrix V — (VJ^i/) 
where 

V^, = i (A4 Ai + a4) = i ({ a4. Ail) (1) 

with A^^ = — (^^) in term of the variables (^^) = 
((7i,Pi, P2) for a two-party system specified by canon- 
ical variables (gi,pi) and (92,^2). We generally define 
(O) — TipO by using the density matrix p. The correla- 
tion matrix is also written in the form 



(2) 



where A and B are 2x2 real symmetric matrices and C 
is a 2 X 2 real matrix. The standard form 



Vo 



a ci 

a C2 

ci 6 

C2 6 



(3) 



is obtained from the general V by applying the 5*^(2, i?)(g) 
Sp{2,R) transformations 3]. 

The separability condition, which is derived from an 
analysis of the non-negativity of the partial transposed 
density matrix p, is written in the matrix notation Q 



with a 2 X 2 simplectic matrix 

J = 



1 

-1 



(4) 



(5) 



or equivalently 



(fAd + fBf -f 2(fCf + g^Ag + h^Bh + 2g'^Ch 



> \d^J9\ + \fJh\ 



(6) 



with / ^ /i all standing for arbitrary real two component 
vectors. 
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When one regards the separabihty condition as a con- 
straint on the range of |ci| and |c2| in the standard form 
Vq in (3), it is written as 



1 



4(a6 - ci){ab - 4) > (a^ + b^) + 2|ciC2| - - 



v/(2a-l)(26-l) > |ci| + |c2| 



(7) 



together with a > 1/2 and b > 1/2. The first algebraic 
relation in (7), which was derived by Simon 0], essen- 
tially corresponds to 



(8) 



and thus it is manifestly invariant under Sp{2, R) (g) 
Sp{2,R). The second condition in (7) is given by the 
weaker conditions derived from (6) with subsidiary con- 
straints g = J^d and h = zLJ^ f, and it is used to ex- 
clude the solutions of the first inequality in (7) which 
allow cf — > oo and 4 oo for fixed a and b. The condi- 
tion a > ^, for example, is given by setting g — j'^d and 
/i = / = in (6). 

The separability condition (7) is explicitly solved, 
namely, the solution of the first inequality which satis- 
fies the second constraint is given by 

c? < ^{[2ab(l + t^)+t]-2^D{a,b,t)}, 

cl < i{[2a6(l + f) +t]- 2^D{a, b, t)} (9) 

with an auxiliary quantity 

D{a, b, t) = a^b^il - t^f + t{a + bt){at + b) (10) 



for 



0<i=|c2|/|ci|<l 



(11) 



where we choose |c2| < |ci| without loss of generality. 
When one defines 

f{cl) = A{ab - cl){ab - t^cj) - {a^ + b^) - 2tcl + J(12) 

for the first inequality in (7), f{c\) assumes a negative 
minimum value at c\ = [2a6(l -f t^) + t]/(4i^) > and 
/(O) = 4(a2 - - i) > 0. Thus /(c?) = has two 

non-negative solutions. Since the second constraint in (7) 
is written as (2a— 1)(26— 1)/(1 + i)^ > c^, we examine 

/((2a-l)(26-l)/(l + t)2) 

-4(a-i)(6-^)H(4-2x) + ix] (13) 

with 0<a; = 4t/(t+l)^ < 1, which is confirmed to be 
negative for < a; < 1 (for the generic case a > \ and 
b > \) and vanish for x — l,i.e.,t = 1. For t = 1 the 



solution in (9) coincides with the second constraint in 
(7) , and thus only the solution in (9) satisfies the second 
constraint in (7). Note that |ci| = |c2| = always defines 
a separable state. 

Another essential ingredient in the proof of Simon Q 
is the P-representation for Gaussian systems. The Gaus- 
sian state is called P-representable if the density matrix 
is written as (we assume (^^) =0 without loss of gener- 
ality) 



P 



d^a / d'^l3P{a,l3)\a,l3){a,P\ 



(14) 



where |a,/3) is the coherent state defined by \a,(3) — 

gaat_i|a|2|Q^ ^^/3bt_i|^|2|Q^ with O = {qi + ipi) / V2 and 

b = {q2 + ip2)/V2. Thus the P-representable states are 
separable. By defining a — (ai + ia2)l\p2 and (i — 
(/?i+i/92)/V2, the weight factor P(a, /3) in (14) is written 
in terms of the correlation matrix V as 



F(a,/3) 
%/detP 



(15) 



1 



47r2 ""^ ' 2 
with the matrix P given by 



cxp{--(ai,a2,/3i,/32)-P(ai,a2,/3i,^2) } 



Y 



1 



-/ > 



(16) 



which defines the condition for the P-representation. See 
also 

The P-representation condition (16) is not invariant 
under 5^(2, i?) ® Sp{2, R) since SS'^ ^ I in general for 
S G Sp{2,R) (g) Sp{2,R). We thus examine the allowed 
range of Sp{2, R) ® Sp{2, R) transformations which pre- 
serve the condition S^^Vo{S'^)^^ — ^I > or equivalently 
Vq — ^SS^ > by starting with the standard form in (3). 
We choose the 5*^(2, R) (g) 5*^(2, R) squeezing matrix S as 



( l/ri 

ri 

l/r2 

V r2 



(17) 



with suitably chosen ri > 1 and r2 > 1. The eigenvalues 
of Vb - \S{ri,r2)S'^ {ri,r2) are given by 



(A 



1 



i)± 



{{a 



2ri 



) + (&- 



2r2 



±^iia-±)-ib-±)r+,cl}, 
(A2)± = ^{(a-iri) + (6-ir2) 



1 



1 



±^Jiia--n)-ib--r2)r+A4}. (18) 

The P-representation exists if (Ai)± > and (A2)± > 0, 
namely, if the following two conditions are simultaneously 



3 



satisfied 



(19) 



1 

2ri 



) + (6 - 277) > and (a - in) + 



together with (a 

{h-\r2) > 0. 

When one regards (9) and (19) as constraints on 
the pair of variables {c\,c^) for given a and &, the P- 
representation condition is more restrictive than the sep- 
arabihty condition; since the P-representation satisfies 
separabHlity, the set of points (cf^c^) allowed by the P- 
representation condition (19) always satisfy the separa- 
bility condition (9). To be precise, we are working on 
the line defined by = cf/cf. See also (30) below. We 
thus expect that these two conditions can coincide only 
for the extremal value of the P-representation condition 
(19) with respect to ri and r2 with fixed t. 

We thus want to prove 



1 



2r2(t,ri) 



) 



^[ia-ln){b-\r,{t,n))] 

1 



= ^{[2ab{l + f) +t]- 2^D{a, b, t)} (20) 

for a suitable 1 < ri < 2a (and 1 < r2 < 2b) for any 
given < t < 1 by regarding r2 as a function of ri and t. 
By this way we establish that the separability condition 
(9) agrees with the P-representation condition (19) with 
a suitable 5*^(2, R) ® 5^(2, R) transformation. 

We first consider the stationary points (or extremals) 
of (Q~ 27r)(^~ 2r2(lri) ) ^nd (a - iri)(6- ir2(t,ri)) in 
(19) with respect to ri with fixed t, namely 



^ 2r2(t,ri)V2 +^'^ 2r,'rldn~ ' 
(&-ir2(t,ri)) + (a-iri)|^-0. 

These two relations combined give rise to 
(ari-1/2) _ (fer2-l/2) 



(a/ri-1/2) (6/r2-l/2)- 



(21) 



(22) 



The relation (22) combined with the first equality in (20) 
gives 



ri = 



n _l_ 112 



^b+ - 
f " ^ 2 

b + 



(23) 



which are symmetric in ri and r2 and are solved as 
1 



ri = 



r2 



at + b 
1 



bt 



{ab{l - <2) + ^D{a,b,t)}, 
{ab{l - t2) -I- ^D{a,b,t)} (24) 



with <t = |c2|/|ci| < 1 and D{a,b,t) defined in (10). 
The squeezing parameters are thus determined. 
One can confirm 



2a > ri > 1, 26 > r2 > 1 



(25) 



for a > ^ and 6 > 5, if one combines the relations (23) 
with 00 > ri/i > 1 and 00 > r2/t > 1 valid for the solu- 
tions in (24). The relation 00 > ri/t > 1, for example, is 
established if one uses the expression 



bt 



n ^ 

t ~ -a6(l - <2) + y/a^b^il - t2)2 + t{a + bt){at + b 

together with t{at + b) < {a + bt) for < < < 1 and the 
triangle inequality. Eq.(24) gives ri = 7'2 = 1 for t = 1, 
and ri — 2a, r2 = 2b for t = 0. 

We emphasize that the condition (25) is required by a 
part of the P-representation condition (a — ^) + {b — 

^) > and (a - ^n) + (b - \r2) > in (19) when 
combined with (22) which implies ^2(1) = 1 and r2(2a) = 
2b if one regards r2(ri) as a function of ri. It should 
be noted that the squeezing parameters themselves are 
essentially determined by the first equality in (20) which 
comes from the P-representation condition. 

We finally evaluate by using ri and r2 in (24) as 



t2 



(a - -ri)(5- -ra) 
= ^{[2a6(l + t2) + t] _ 2^D{a, b, t)} (27) 

which is a remarkable identity. This relation establishes 
(20), namely, the fact that the boundaries of the condi- 
tions for separability and P-representation coincide for 
any < t = |c2|/|ci| < 1. It is significant that the 
squeezing parameters in (24) depend only on the ratio 
t — \c2\l\c1\ and not on |ci| and \c2\ separately unlike 
the case in [(^; this is because we are determining the 
bound to |ci| for given t to be consistent with (9). All the 
states parameterized by |ci|, which satisfy (9) for given 
a, b and t, automatically satisfy the P-representation 
with the parameters in (24). 

Our explicit construction proves that the P- 
representation condition (19) with the squeezing param- 
eters in (24), which satisfy 1 < ri < 2a and 1 < ^2 < 26, 
is equivalent to the separability condition (9) for any 
< t = |c2|/|ci| < 1, and thus the separability condi- 
tion (9) is a necessary and sufficient separability crite- 
rion for two-party Gaussian systems. To our knowledge, 
our study is the first quantitative treatment of the basic 
algebraic condition (7) which is considered to be funda- 
mental The existence of explicit analytic solutions 
of ri and r2 in (24) for this basic problem is interesting, 
and the triplet 



(Vb,ri,r2) 



(28) 



where Vb is the standard form in (3), characterize the 
general covariance matrix (2) of either separable or insep- 
arable Gaussian states. These explicit solutions are con- 
venient for quantitative treatment, and thus they should 
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be useful for the quantitative theoretical and experimen- 
tal analyses of entanglement in practical applications. 



III. IMPLICATIONS OF ANALYTIC 
SOLUTIONS 

We would like to compare the present analysis with 
some of the past related works. The separability con- 
dition (4) is based on the partial transpose opera- 
tion of the density matrix in the manner of Peres and 
Horodeckidi]. The analysis of the negativity of the par- 
tial transposed density matrix has been further extended 
by Shchukin and Vogel 0, • We here note a complemen- 
tary property that the second moment for the separable 
density matrix p = Y^nPn^^^n '^ith ipn = 4>niqi)'Pn{q2) 
gives rise to a generalization of (6) 

(fAd + fBf + 2(fCf + g^Ag + h^Bh + 2g^Ch 

> (fAd + fBf -I- 2d^Cf + g^Ag + h^Bh + 2g^Ch 
+ \d^Jg\ + \fJh\ (29) 

by an analysis similar to the derivation of the Kennard's 
uncertainty relation without referring to the partial trans- 
pose operation. Here V in terms of A, B and C is de- 
fined analogously to (2) with the elements of the ma- 
trix V = (V^,) defined by V^, = EfcPfe(A4)fc(Ai)fc 
and A^^ = — (Cm)- Note that one may choose 
(Im) = I]fePfe(4)fe = 0' but the average (|^)fc for each 
component state does not vanish in general. A and B 
are 2x2 real symmetric matrices and C is a 2 x 2 real 
matrix. Both of V in (2) and V are non-negative, and 
thus the condition (6) is necessary but not sufhcient for 
separability in general. For the P-representation of the 
Gaussian state in (15) one can identify P^^ = by us- 
ing a special property of the coherent state, which spoils 
Sp{2, R) (X) Sp{2, R) invariance due to normal ordering, 
and at the boundary of the P-representation condition 
where two eigenvalues of vanish, the condition (29) 
can be equivalent to (6). [§]. 

One can also directly derive the condition (6) from the 
P-representation condition (16) which implies 

d^Ad + fBJ + 2d^Cf -\{d^d + ff) > (30) 

for any d and /. When one adds (30) to a relation ob- 
tained from (30) by replacing d and f hy g and h, respec- 
tively, one recovers the separability condition (6) 

d^Ad + fBf + 2d^Cf + g^Ag + Bh + 2g'^Ch 

> ^{d^d + ff) + i(.g^5 + h'h) > \d^Jg\ + \fjh\, 

where we used the relation such as ^{d^d + g^g) > 

^idTd){gTjTjg) > \dTjg\. 

What we have proved in the present paper (and also 
in 0]) is that (6) also implies (16) with the help of the 



squeezing operation. In this proof, the second weaker 
condition in (7) played a crucial role, namely, the first 
relation in (7) does not encode the entire information of 
the separability condition (4) or (6). The importance 
of the second condition in (7), which appears to be over- 
looked in the past related works (see, for example, eq.(19) 
in [3!], Theorem V.2 in Q and eq.(23) in 0), was first rec- 
ognized clearly in the present explicit analytic treatment 
of the separability condition (6). 

As for the second weaker condition in (7) , which is not 
Sp{2, R) Sp{2, R) invariant, one can confirm that this 
weaker condition 



C2 



:(31) 



an + — - l][6r2 + — - 1] > V^lcil , , 



written for the second moment, which is obtained from 
the standard form Vq in (3) by a squeezing operation 
in (17), corresponds to the separability condition in 
in fact, if one imposes the condition (22), the left-hand 
side of (31) becomes 



a nr. b , 

[an + l][br2 + 1] 

ri r2 



(32) 



= v/(ari - 1/2) (6r2 - 1/2) + ^{a/n - 1/2) (6/^2 - 1/2) 

and one recovers eq. (16) in [6^ when converted into their 
notation. 

From our analysis of (13), it is obvious that some 
parameter region allowed by the weaker condition (31) 
with ri = r2 = 1 does not satisfy the P-representation 
condition. But if one combines this weaker relation 
with suitable squeezing, a simple direct proof of the P- 
representation is obtained. If one sets |c2| = t\ci\ in (31) 
with the condition (22), one obtains 



v/(ari - 1/2) (6r2 - 1/2) + ^{a/n - 1/2) (6/^2 - 1/2) 



]\ci\ 



(33) 



If one uses the analytic formulas of squeezing parame- 
ters given in (24), one can confirm that the relation (33) 
when regarded as a bound to |ci | becomes identical to the 
Sp{2, R) (g) 5*^(2, R) invariant bound (9) and also to the 
boundary of the P-representation condition (19). Here 
we use (20) and its square root, namely. 



V(ari^l/2)(&r2-T72) 

^ V(a/ri-l/2)(6/r2-l/2) 

= l{[2a5(l + i2)+t]_2/DKM)}'/'. (34) 

Since (19) implies (33), one concludes that the weaker 
separability condition combined with our explicit squeez- 
ing parameters provides the necessary and sufficient con- 
dition for the P-representation. In a symbolic notation 
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we have 



eq.(33) D cq.(9) D eq.(19), 



(35) 



which means, for example, any standard form of the co- 
variance matrix Vq which satisfies (19) automaticahy sat- 
isfies (9) (a stronger condition means a smaller set of Vq). 
But our analysis above shows that (33) and (19), both of 
which depend on squeezing parameters, coincide with (9) 
for our explicit solutions of squeezing parameters. This 
illustrates the power of our explicit analytic formulas, 
and this simple proof of the P-representation gives an- 
other explicit proof of the necessary and sufficient con- 
dition for separability of two-party Gaussian systems on 
the basis of (33). 

We next briefly discuss the quantitative difference be- 
tween our scheme and the scheme in 6] . The authors in 
look for the solution /(r^) — where 



fin) = [Vn^\ci\ - Vian - l/2)(&r2 - 1/2)] (36) 
- [|C2|/V^- VNn - l/2){b/r, - 1/2)] 

and the condition (22) when written in our notation. 
Namely, they look for the solution of 

y/{an - l/2)(fer2 - 1/2) - V^|ci| 

= V(a/n - l/2)(6/r2 - 1/2) - \c2\/Vn^ (37) 

together with (22) and (25), though the condition (25) is 
not explicitly stated in 6]. In contrast, in our scheme we 
solve 

v/(ari-l/2)(fer2-l/2) 



^/nr2\ci\ 
v/(a/ri-l/2)(&/r2-l/2) 

i\c2\/y/n^) 



(38) 



together with (22) and (25), as is seen from the first 
equality of (34). Since (38) implies 

V(ari-l/2)(6r2-l/2) - V^l^i I 



y/rir2\ci\ 



V(a/ri-l/2)(6/r2-l/2)- \c2\/V^ 



{\C2\/Vn^) 

the common solutions of (37) and (38) exist only for 



, (39) 



v'(ari-l/2)(for2-l/2) - Vn^\ci \ (40) 
= \/(a/ri-l/2)(6/r2-l/2) - \c2\/Vn^ = 



or else for ^/rlrj|ci| = |c2|/^/nTJ, namely, 

\C2\ 



rir2 



= t. 



ci 



(41) 



The relation (40) together with our explicit analytic so- 
lutions of ri and r2 shows that the values of |ci| and \c2\ 



given by (40) correspond to the largest allowed values of 
|ci| and |c2| in the separability condition (9) if one re- 
calls (34). As for (41), we can choose < t < 1 without 
loss of generality as in (11) and the squeezing parameters 
are bounded as in (25) by the P-representation condition 
when combined with (22). Thus the condition (41) is sat- 
isfied only at t = 1 and ri = r2 = 1. The two conditions 
(37) and (38) thus coincide only for such a measure-zero 
subset of the standard form of separable covariance ma- 
trices {Vo}. This quantitative difference between the two 
schemes is interesting. 

As we explained already, the squeezing parameters in 
our scheme specify the boundary of the P-representation 
condition for any given t, which agrees with the separa- 
bility condition (9). On the other hand, the scheme in ^ 
specifies the squeezing parameters for each given covari- 
ance matrix separately such that the P-representation 
condition is ensured, and thus the specification of squeez- 
ing parameters is more local. Our analytic solutions de- 
pend only on the ratio t = |c2|/|ci|, while squeezing pa- 
rameters in (37) generally depend on |ci| and |c2| sepa- 
rately. Nevertheless, we here show that our explicit ana- 
lytic solutions allow a concrete proof of the separability 
criterion in the scheme of Q . For this purpose, we extend 
(35) to a symbolic notation 

eq.(33) D eq.(9) D eq.(19) D {eq.(33) n eq.(37)} (42) 

where the last relation means that any Vq which satisfies 
(37) in addition to (33) automatically satisfies the P- 
representation condition (19), as is explicitly confirmed 
^6']. Thus {eq.(33) fl eq.(37)} provides a sufficient condi- 
tion for the P-representation, but the converse is not ob- 
vious. The P-representation condition (19) implies (33) 
as is shown in but it is not obvious if (19) implies (37). 
To be more precise, it is not obvious if (9) impfies (37). 
In fact, the authors in Q choose a rather general class 
of Vo and prove the solution of f{r*) = in the domain 
1 < r* < oo which, however, does not satisfy a part of the 
P-representation condition (25), namely, 1 < r* < 2a. 

We now show that (9) implies (37) by using our ana- 
lytic solutions. We first write (36) in the form 



/(r i,|ci|)^(V?I?^-WF I?^)|ci| ^ 

-^(ari - 1/2) (6r2 - 1/2) + ^{a/n - l/2)(6/r2 - 1/2) 

together with (22). We then have /(ri = l,|ci|) > 
for < t < 1 (the case t = 1 is trivial since then ri — 
r2 = 1 is the solution). From the analysis of (40), we 
have /(ri, |ci|) = for the largest allowed |ci| in (9) and 
for our analytic solutions. Since y^rlrj > 1, the function 
/(ri, |ci|) is a linear increasing function of |ci| when we 
fix squeezing parameters at our analytic solutions. Here 
it is crucial that our analytic solutions (24) depend only 
on a, b and t. We thus conclude 



/(^i,|ci|) <0 



(44) 



for all values of |ci| in (9) when the squeezing parameters 
are fixed at our analytic solutions. Since our analytic so- 
lutions satisfy 1 < ri < 2a as in (25), the relation (44) 
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combined with /(ri = 1, |ci|) > shows that f{r^) = 
has a solution in the interval 1 < r* < 2a for all values 
of |ci| in (9). Namely, (9) imphes (37). This completes 
the proof of the necessary and sufhcient separability con- 
dition in the scheme of [(|. 

In the above analysis, we used the solution of (22) 



[^{1 - Xf + IGb^X + (1 - X)] 

with X{ri) — {2a/ri — l)/(2ari — 1) which assumes 
X{1) = 1 and X{2a) = 0, and thus r2(l) = 1 and 
r2(2a) = 2b. 

IV. DISCUSSION 

We found the explicit analytic formulas of squeez- 
ing parameters which establish the equivalence of the 
Sp(2, R) (X) 5*^(2, R) invariant separability condition with 



the P-representation condition. These explicit analytic 
formulas give rise to not only the concrete proofs of the 
separability criterions for two-party Gaussian systems 
formulated in J| and ^ but also a new simple proof 
of the separability criterion as is described in connection 
with (35). Our analytic formulas thus provide a unified 
basis for the analysis of separability in continuous vari- 
able two-party Gaussian systems. 

We analyzed the two-party system with one freedom 
for each party. The system we analyzed may be more 
properly called a two-mode system since the analysis of 
the two-party system with more than one freedom in each 
party is more involved [lol |. As for the non- negativity 
of the partially transposed density matrix, Vidal and 
Werner analyzed the separability of Gaussian states 
by using logarithmic negativity, which is essentially the 
same as the actual analysis of Simon ^ and leads to (7) 
and (16), but no analytic formulas of squeezing parame- 
ters are given. See also the related references [ll] - [l9| . 
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